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ABSTRACT 

A recent analysis of the relationship between optical and X-ray properties for the 
Einstein sample of early-type galaxies has revealed that SO galaxies have lower mean 
X-ray luminosity Lx per unit optical luminosity Lb than do ellipticals. In the same 
analysis, significant correlations are found between the X-ray properties and the axial 
ratios, such that the roundest systems have the highest Lx/Lb', this trend holds also for 
either Es or SOs alone (Eskridge et al. 1995a, b). The systematic X-ray underluminosity 
of SOs with respect to Es could be explained with a higher heat input or with a 
lower gravitational energy (all per unit gas mass), at fixed Lb- The gravitational 
energy could be lower because of their higher rotation, which decreases the effective 
potential, or their different mass distribution. These possibilities are examined here, 
by considering their role in the global energy budget of the hot gas flows in early- 
type galaxies. The effect of the flattening of the mass distribution is investigated 
with galaxy models described by the Miyamoto-Nagai potential-density pair, to which 
a dark matter halo of various shapes is added. For these two-component models 
the analytical expressions of the gravitational energy, and the stellar kinetic energy 
associated with various relative amounts of random motions and rotational streaming, 
are given. It is found that rotation cannot produce a change in the flow phase of the 
hot gas, independently of the galaxy shape and the presence of dark matter. The 
effect of flattening instead can be substantial in reducing the binding energy of the 
hot gas. Thus SOs and possibly non spherical Es are less able to retain a significant 
halo of hot gas than rounder Es of the same Lb ■ 

Key words: galaxies: elliptical and lenticular, cD - galaxies: ISM - galaxies: kine- 
matics and dynamics - galaxies: structure - X-rays: galaxies 



1 INTRODUCTION 

X-ray observations, beginning with the Einstein Observa- 
tory (Giacconi et al. 1979), have demonstrated that normal 
early-type galaxies are X-ray emitters, with 0.2-4 keV lumi- 
nosities ranging from ~ 10^'' to ~ 10^^ erg s~^ (Fabbiano 
1989; Fabbiano, Kim, &; Trinchieri 1992). The X-ray lu- 
minosity Lx is found to correlate with the blue luminosity 
Lb {Lx <x L^^^^'^), but there is a large scatter of ap- 
proximately two orders of magnitude in Lx at any fixed 
Lb > 3 X 10^'' i0. The observed X-ray spectra of the 
brightest objects are consistent with thermal emission from 
hot, optically thin gas (Canizares, Fabbiano &; Trinchieri 
1987). These hot gas coronae are accumulated during the 
galaxy lifetime from the stellar mass loss, and are heated 
to X-ray temperatures by the thermalization of the stellar 
random motions, and by the type la supernova explosions 



(hereafter SNla). The scatter in the Lx — Lb diagram has 
been explained in terms of environmental differences (i.e., 
varying degrees of ram pressure stripping due to the inter- 
action with the intracluster or intragroup medium. White &; 
Sarazin 1991), or in terms of different dynamical phases for 
the hot gas flows, ranging from winds to subsonic outflows 
to inflows (Ciotti et al. 1991, hereafter CDPR). 

Recently Eskridge, Fabbiano &; Kim 1995a,b conduct- 
ed a multivariate statistical analysis of data measuring the 
optical and X-ray properties of the Einstein sample of early- 
type galaxies; this includes 72 ellipticals and 74 SOs, and is 
the largest X-ray selected sample of such galaxies present- 
ly available. Eskridge et al. 1995a show that on average SO 
galaxies have lower X-ray luminosity at any fixed Lb than 
do ellipticals (see their Fig. 8a,b,c); a comparison of the 
distribution functions for the two morphological subsamples 
indicates that the SOs have lower mean Lx (at the 2.8a level) 
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and Lx/Ljj (at the 3.5a level) than do the Es. Morevover 
Eskridge et al. 1995b find that galaxies with axial ratio close 
to unity span the full range of ixi while those with large ax- 
ial ratio all have ix^lO^^ ^^g s~^. The relationship defined 
by ix/^B is stronger than that defined by ixi it is in 
the sense that those systems with the largest ix/^B ^^e the 
roundest; this means that at any fixed ig the X-ray bright- 
est galaxies are also the roundest. This correlation holds for 
both morphological subsets of Es and SOs. Thus SOs and 
non spherical Es seem to be less able to retain hot gaseous 
halos than are rounder systems of the same ig . This could 
be due to a lower gravitational energy, because of their high- 
er rotation, which decreases the effective potential, or their 
different mass distribution (see also Pellegrini 1994). 

The purpose of this paper is to investigate whether a 
flat and partially rotationally supported galaxy is expected 
to host a different gas flow phase with respect to a spherical 
pressure supported galaxy of the same ig . This is accom- 
plished with the study of the energetic balance of the hot 
gas flows. Axisymmetric two-component galaxy models are 
built, where the stellar and dark mass distributions have 
the Miyamoto-Nagai shape; the shape of the dark matter 
halo is varied from flat to spherical. The decoupling be- 
tween ordered and disordered azimuthal motions is assumed 
to satisfy the Satoh's ^-decomposition, or a "maximum ro- 
tation" decomposition. The analytical expressions for the 
virial terms of these two-component models, together with 
the projected kinematical properties for the one-component 
model, are given here for the first time. 

In Section 2 the role of flattening and rotation in the 
energetic balance of a galaxy of a general shape is quali- 
tatively studied. In Section 3 the two-component galaxy 
models are introduced; in Section 4 and 5 these are used 
to derive quantitative estimates of the roles of rotation and 
flattening. In Section 6 other possible explanations for the 
X-ray underluminosity of flat systems are discussed, and the 
main conclusions are summarized. 



2 THE ENERGY BUDGET OF THE HOT GAS 

For pressure supported spherical galaxy models, a reliable 
estimator of the nature of the flow phase of the hot gas can 
be derived from its global energy balance (CDPR). This es- 
timator is the ratio % between the power required to steadily 
extract the stellar mass loss from the galaxy potential well, 
and the power supplied to the gas by the thermalization 
of the stellar random motions and by the supernova heat- 
ing. Spherically symmetric hydrodynamical simulations, in 
which radiation cooling is taken into account, show that in 
the scenario of subsequent wind, outflow, and inflow (here- 
after WOI), X describes in an accurate way the nature of 
the flow. The simulations show the correspondance between 
X < 1 and the presence of a wind phase (i.e., of a low £x)i 
and the onset of a high ix soon after x becomes > 1 (CD- 
PR). They also show that small variations (within 10%) of x 
when it is close to unity describe large flow regime variation- 
s, from winds to inflows. The situation is obviously different 
in other scenarios where x is not < 1 from the beginning of 
the flow history, i.e., if the initial flow phase is not a wind 
(and so a major fraction of the heating goes into radiation). 
The use of x is extended here to the case of galaxies of a 
general shape and internal dynamics. 



2.1 The generalized parameter x 

For the needs of the following discussion, the classical scalar 
virial theorem for a stellar distribution of spatial density p* 
interacting with a dark matter potential $h is recalled: 

2T + n = |W^|, (1) 

where H = J p^,TT[a^)dV is twice the kinetic energy 
associated with the stellar random motions, Tr((j^) is 
the trace of the velocity dispersion tensor of the stars; 
T = 0.5 J p^,^I[v[Vj)dV is the kinetic energy associated 
with the stellar ordered motions; \W\ = |VFh,| + |Wh| = 
— J Ph.(0.5$h. — X]idiYi/ dx]i)dV , where is the interac- 

tion energy of the stars with the dark halo. 

The terms entering x now defined. The power Lg 
required to steadily extract from the galaxy the stellar mass 
loss is 



I 



p4^* + ^i,)dV=\U\ = 2\W,\ + \Ui,\, 



(2) 



where a(i) = M^/M^ is the stellar mass loss rate Af* per 
unit stellar mass, Af* being the total stellar mass (see CD- 
PR). The thermalization of the stellar random motions heats 
the gas at a power given by 



4^= 1 [ p,TT{a')dV=^. 
a{t) 2j''^> 2 



(3) 



Indicating the SNIa heating per unit time with £sNi the x 
corresponding to a non rotating, totally velocity dispersion 
supported galaxy is: 

xn = ^^= . J^L_, . (4) 



isN/a + 0.5|VF| 



Note that formally < xn < 2|?7|/|W^|, and that xn in- 
creases going from a very high to a null SNIa heating. 

For a fixed mass distribution (i.e., a fixed |W^|), the 
amount of rotational streaming can formally vary from zero 
to |VF|/2; the power irot related to this fraction of the total 
kinetic energy is given by 



p^,T^(v{v^)dV = T. 



(5) 



-^rot 
a(i) 

How does irot enter eq. (4)? There are two extreme views. 
In the first one, the ordered motion is considered as a change 
in the effective potential experienced by the gas: in this 
case irot is to be subtracted to the numerator of (4). In 
the opposite view all the kinetic energy of the gas, also the 
ordered one, is eventually thermalized: in this case irot is 
to be added to the denominator of eq. (4). We parameterize 
the real behavior - which lies between the two extreme cases 
- by defining 

ig - T-trot \U\ - jT 



X ■ 



isN + io- + (1 - 7)ir 



isN/a + 0.5|VF| -7T 



(6) 



with < 7 < 1. Note that if 7 = then x = Xlli i-^-i if the 
thermalization of irot is complete, x coincides with xn even 
when rotation is present. When instead 7 = 1, i.e., in the 
case of no thermalization, the effect of rotation is maximum. 

In the degenerate and trivial case when rotation, dark 
halo, and SNIa's are absent, x = 4 from eqs. (6) and (2); the 
same applies when rotation is fully thermalized, and again 
the dark halo and the SNIa's are absent. In all cases when 
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the supernova heating is negligible, x > 1 an inflow is 
present. 



2.2 The role of flattening and rotation 



We consider now two galaxies of the same luminosity ig, 
i.e., with the same SNIa heating, and the same stellar mass. 
One of these is the reference galaxy, a completely velocity 
dispersion supported system; the other galaxy is the test 
object, with which the effects of flattening and rotation are 
investigated. For this study it is useful to consider the ratio 
between % of the test galaxy and xn of the reference one, 
whose properties are characterized by the subscript "0": 



X 

xn 



jT 

\U\ 



\u\ 



+ xn 



\W\ 



\Wo\ 



2\U\ 



\U\ 



(7) 



A general variation in the structure and internal dynamics 
of the reference galaxy can be obtained through two inde- 
pendent steps: the first one changes the structure from the 
reference case to the test one, and the second step adds the 
rotation. The first step implies turning Uq and Wq into 
U and W, e.g., changing the stellar mass distribution and 
dark matter amount and/or distribution. From eq. (7) with 
T = a dimensionless parameter is found 



Q 



\Uo\ 



|[7|_0.5(|W^|-|W^o|) 



(8) 



whose value determines the resulting %. When, by changing 
the galaxy structure, Q becomes < xjji then % > 1; when 
Q = xn, then X = 1; and when Q > xni then x < 1- 

The second step adds rotation to the modified galaxy. 
Using again eq. (7) with \W\ = \Wo\ and |?7| = |?7o|, and 
with Xn the resulting x of the first step, it results that: 
1) if xn = li also X = li for any jT; 2) if xn > li also 
X > 1, and an increase in jT produces an increase in xi 
3) if Xn < 1 then x < li and an increase in jT produces 
a decrease in x- In other words, the effect of rotation is 
never that of changing a x > 1 into a x < li and vicever- 
sa, independently of the unknown value of 7. A qualitative 
explanation of this behavior can be found by inspecting e- 
q. (6): when the gas in absence of rotation is bound, the 
introduction of rotation causes the subtraction of the same 
quantity from the numerator and the denominator, and this 
corresponds to a larger percentual decrease of the heating 
than of the binding energy. So, an increase in x is produced. 
The reverse is true when the gas is unbound in absence of 
rotation. As a consequence, in the WOI scenario, rotation 
cannot produce a change in the flow phase; instead, it acts 
in the sense of making even more stable the existing flow 
phase, for any fixed galaxy structure. It would be inter- 
esting, for the purpose of the present study, to know under 
which conditions and how frequently, at ig fixed, the flat- 
tening produces Q > xni where xn refers to a spherical 
galaxy and is > 1; in this case in fact x would be < 1 for 
the flattened object, and this could host an outflow. If also 
some rotation is added, a further decrease in x is produced 
(if 7 > 0), and the flow is even more likely to be an outflow 
or a wind. In this way the average X-ray underluminosity 



Figure 1. The variation of x produced by structural changes. 
Top panel: no dark matter. In the middle and lower panels, for 
the same xu values, solid lines correspond to s = 0, dotted lines 
to s = 1, and dashed lines to s = 5. In the middle panel the 
dark matter keeps its spherical shape while the stellar mass is 
flattened; in the lower panel sj^ = 0.5s. 

of flat systems at fixed ig would be explained. What is 
the efficiency of the flattening in producing a variation in Xi 
for realistic galaxy models? What is the importance of the 
further variation that can be produced by rotation? These 
points will be investigated in the case of SOs in Sections 4 
and 5, using the galaxy models described below. 



3 THE MODELS 

The stellar density distribution of an SO galaxy is mod- 
eled by adopting the potential-density pair introduced by 
Miyamoto &; Nagai (1975), as given in Binney &; Tremaine 
(1987, p.44; hereafter BT87): 

'•^^ gJl^ + (g + 3C)(a + C)^ 



p* 



47r J ^3[JJ2 +(a + ^)2]5/2 



Vii2 + (a + C)2' 



(9) 



(10) 



where ^ = + b'^ , and [R,ip,z) are the cylindrical coor- 
dinates. For a = the stellar density is a Plummer (1911) 
sphere with effective radius b; for 6 = eq. (9) describes the 
Kuzmin (1956) disk. A constant stellar mass-to-light ratio 
is assumed throughout the system. The Miyamoto-Nagai 
(hereafter MN) model produces a realistic gross structure, 
similar to that observed for SOs; particularly, the light dis- 
tribution resembles that of disk galaxies when a/b = 5 (e.g., 
BT87). In order to study the effect of a dark matter halo, 
another Miyamoto-Nagai density distribution characterized 
by Afhi '^hi and fc^i is added to the stellar profile. 

The quantities entering x are determined using the 
Jeans equations. Following Binney, Davies, &; lUingworth 
(1990) we implicitly adopt a distribution function that de- 
pends only on the isolating integrals E and Lz- In this case: 
1) the terms v^Vz, v^Vi^, and Vi^Vz are zero (and so the veloc- 
ity dispersion tensor is aligned with the coordinate system); 
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2) the radial and axial velocity dispersions are equal, i.e., 
(jR = (Jz = (J; 3) the only nonzero streaming motion is in 
the azimuthal direction {f). The resulting Jeans equations 
are: 

,2 



dz 
and 

OR 



-p* - 



+ p* 



dz ' 



R 



OR 



(11) 



(12) 



where $tot = $* + $h- 

Since in eq. (12) only v'^ appears, the decomposition 
of the azimuthal velocity can be chosen freely. Being the 
stellar density flat, a natural choice is the ^-decomposition 
introduced by Satoh (1980): 



= k {v^ — )i 



and then 



(13) 



(14) 



with < i; < 1. In the case i; = all the flattening is due 
to the azimuthal velocity dispersion, i.e., no net rotation is 
present; < i; < 1 is a sufficient condition in order to have 
a positive value of cr^ everywhere, and k = 1 corresponds 
to the isotropic rotator. This solution does not describe 
the maximum rotation formally allowed for a given galaxy 
structure; this is obtained by relaxing the condition that k is 
a constant, and requiring that cr^ = over the whole space: 

-• (15) 



*:max(-R,^) = 1 + = 



While certainly such a solution is not physically possible for 
a realistic structure, nevertheless it is used as a limiting case 
in the investigation of the maximum effect of rotation. 



4 MODELS WITH NO DARK HALO 

In the case of no dark halo, the Jeans equations can be solved 
analytically (see Nagai &; Miyamoto, 1976). The velocity 
dispersion is recovered from (11): 



p*cr 



2 _ GM^b' (g + C) 

Stt C2[-R2+(a + C)2p' 



and after some easy algebra one obtains from (12): 

P*{V^ - (J ) ^ 



R' 



(16) 



(17) 



(18) 



A-K C^[R-i + {a + CYf 
Using eqs. (15), (16), and (17) we can recover 

^max(^i^) = 1 H 2aR? 

The analytical expressions for the various terms entering 
the virial theorem have been derived, both for the constant 
k case and the maximum rotation case (see Appendix B). 

4.1 The effect of flattening and rotation 

We compute here the % variation produced by the flattening 
and the rotation for MN models; the galaxy total luminosity 



is kept constant. As outlined in Section 2.2, the effect of 
these two processes can be considered independently, and 
we start examining the effect of flattening onto a spherical 
reference galaxy. Two structural parameters describe the 
process of flattening: s = a/b, and r = b/bg, where bg is 
the effective radius of the reference galaxy, s measures the 
intrinsic shape of the flattened galaxy, and r measures its 
size. 

The parameters s and r enter eq. (7) only through Q, 
and the effect of flattening is given by 

4-Q 



(19) 



X ^ 

xn 3Q + xn(i-(3)' 

Note that, in case of no dark matter, Q = 4/(1 + 3ig/ig^o): 
so Q is directly related to the variation in the potential ener- 
gy of the hot gas. Figure 1 (upper panel) shows % vs. Q, for 
three chosen values of xn- A change in the galaxy structure 
is able to produce large variations in %, of the order of those 
needed for a change in the flow phase. More precisely, since 
for spherical galaxy models xn hardly exceeds 1.5 (CDPR), 
at most Q of this order is required to have % < 1. 

The relationship between Q and {s,r) is shown in Fig. 
2, where the contour levels of constant values for Q[s,r) are 
plotted. The interesting region in the {s,r) plane is that in 
which Q[s,r) > 1: only in this case a system with xn > 1 
might become characterized by % < 1 after a flattening, 
if this produces Q > xn- In other words, given a spherical 
galaxy with xn > Ij all (s, r) pairs lying in the upper part of 
the plane which has as lower bound the curve Q{s, r) = xn 
describe flattened galaxies with % < 1. It is interesting 
to note that an increase in Q (i.e., a decrease in ig) for 
roughly spherical galaxies is more easily obtained through 
their flattening than through an increase in their size, while 
the reverse is true for disky objects. 

In order to find the location of real galaxies in the (s, r) 
plane, in Fig. 2 two families of lines are shown: the first 
one describes the ratio of the projected central velocity dis- 
persion (jp of the fiattened system seen face-on to the same 
quantity for the corresponding spherical reference galaxy [e- 
q. (A15)]; the second family describes the same ratio using 
the edge-on velocity dispersions [the choice of the decompo- 
sition for the azimuthal component of the velocity dispersion 
cannot affect the central value of the projected velocity dis- 
persion, as can be checked with eqs. (AS), (A9) and (A12)]. 
From Fig. 2 two important results are derived: 1) reason- 
able values for r and the (jp ratio (i.e., near unity), and for s 
(<1), correspond to Q values that can produce % variations 
of the order of those required to have a change in the fiow 
phase; 2) the region where Q > 1 is also the region where 
the fiattened systems have a lower central velocity dispersion 
than the spherical systems of the same luminosity. 

A qualitative explanation for the behavior of Q and 
the central velocity dispersion can be given. The pro- 
jected velocity dispersion is proportional to the average 
of the potential over the stellar density, and, as a con- 
sequence, at the zeroth-order the central (jp contains in- 
formation on the v alue of t he centr al poten tial. It fol- 
lows that (Jn oc ^-$H.(0) oc l/^7-(l + s). Similarly, 
ig OC — Ph.(0)$h.(0) oc (3 + s)/7'^(l + s)^, and so increasing s 
and r lowers (jp and much more ig. This qualitatively ex- 
plains why Q variations that are very effective in producing 
a change in the fiow phase correspond to small changes in 
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Figure 2. Contour levels of constant Q values (solid lines), in 
the case of no dark matter. Dashed lines show the loci where the 
ratio between the edge-on (Tp of the flattened and the spherical 
system is 1, 0.75, 0.5; dotted lines describe the same ratios for 
the face-on projection. 



(jp. Eskridge et al. 1995b find a significant trend of the axial 
ratios with (jp , in the sense that galaxies with the largest (jp 
are the roundest, for constant ig. So, a scenario in which % 
is preferentially < 1 in flat systems can explain their system- 
atic X-ray underluminosity, and finds support in the optical 
data. 

As outlined at the end of Section 2.1, the previous re- 
sults apply also when rotation is fully thermalized. The 
efficiency of rotation in producing a variation in % is investi- 
gated now for the limiting case of no thermalization (7 = 1), 
i.e., when its effect is maximum. Two cases are considered: 
the isotropic rotator, T = T;s [i; = 1 in (B9)], and the max- 
imum rotation solution, T = Tmax, given by (Bll). Figure 
3 (upper panel) shows the variation in % produced by ro- 
tation, for a fixed galaxy structure characterized by various 
degrees of fiattening. This variation is always very small: 
for example if, after a fiattening with s = 5, xn = O-^i then 
the additional reduction in the case of the maximum rota- 
tion solution is 3%. So, not only rotation cannot change the 
gas fiow regime, but its effect is almost negligible anyway. 

5 MODELS WITH DARK HALO 

The analysis of the previous section is extended here to the 
two-component MN models. For a given stellar structure, 
the parameters defining the dark halo are: 72. = Af^/Af*, 
6hi and sj, = a^/fc, and so the most general change in the 
total mass distribution involves five parameters. Note that 
the freedom in the choice of aj, permits to change the fiat- 
tening of the dark mass independently of that of the stellar 
mass. In order to reduce the number of free parameters, it 
is assumed that 6^ = b, which allows to perform the cal- 
culations analytically; moreover, 72. is kept constant during 
the process of fiattening, i.e., round and disky objects of the 
same ig have the same amount of dark matter. Finally, it 
is assumed that the dark matter dominates the dynamics, 
i.e., in eqs. (7) and (8) only the contribution of the dark 
matter potential is considered. So, three parameter values 
are to be specified: r, s, sj,- 

In Appendix B the analytical formulae for the terms 
entering the virial theorem and % are given, for the Satoh's 
and the maximum rotation decomposition of the azimuthal 
velocity. 



Figure 3. The variation of x produced by rotation; x and xu 
refer to a galaxy with the same structure. The various lines cor- 
respond to various degrees of the flattening for the stellar mass 
(s = 0, 1, 5), and larger variations correspond to larger s. In the 
upper panel there is no dark matter, in the middle the dark halo 
is spherical, in the bottom sj^ = 0.5s. 

5.1 The effect of flattening and rotation 

As previously done in the case of no dark matter, we start 
studying the trend of % with the fiattening: 

^ = i^^^^ ^ (20) 

(4-xn)Q + 2(xn-Q)W^h(c,^)/t/h(c) 

where c = [s + Sh)/2, and the functions and fJ^ are 
given in (B3) and (B6). In Fig. 1 we plot % versus Q for 
three fixed values of xn and s. Concerning the shape of 
the dark halo, two cases are investigated: in the first one 
the dark matter distribution is spherical for any fiattening 
of the visible matter (sj, = 0), in the other the process of 
fiattening affects also the halo mass distribution (sj, = 0.5s). 
The case = s coincides with the case of a one-component 
MN model. Figure 1 shows that the presence of the dark 
matter does not change significantly the trend shown pre- 
viously: the same values of Q produce essentially the same 
variations in % that were produced in the case of no dark 
matter. The relationship between Q and {s,r) is very sim- 
ilar to that obtained in the no dark matter case (see Fig. 
4). 

Again we overplot in the {s,r) plane the contours of 
constant value for the ratio of the face-on (Jp's in the fiat 
and in the spherical reference galaxy. When sj, = 0.5s, con- 
clusions very similar to those derived without dark matter 
are valid. As before, the Q contours have roughly the same 
shape of the (jp contours, but the relative position of the 
two has changed slightly: now the (jp contours are shifted 
upwards with respect to the Q contours. As a result, for 
example, the line of equal value for (jp in the spherical and 
in the fiat system lies above the Q = 1 line, and so even fiat 
systems with (jp comparable to that of the spherical galaxy 
can have % < 1. 
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Figure 4. Contour levels of constant Q values, in the case of dominant dark matter (solid lines). In each panel the dotted lines show 
the loci where the ratio between the face-on (Tp of the flattened and the spherical system is 1.25, 1, 0.75. 



(jp shows instead a clearly different behavior from the 
case without dark matter, when the dark halo is kept spher- 
ical (Fig. 4, left panel). Note in particular that high Q 
values can be reached while keeping (jp similar in the spher- 
ical and in the test galaxy. A qualitative explanation of the 
behavior of Q and (jp follows from considerations similar 
to those of Secti on 4.1. When th e dark matter dominates, 
(Jp oc 1/ ^(l + the potential energy 

of the gas ig oc -pH.(0)$h(0) oc (3 + s)/7-''(l + s)^{l + Sh)- 
If Sh = 0.5s, increasing s lowers (jp and ig in the same way 
as in the no dark matter case. On the contrary if s increases 
when and r are fixed, the (jp ratio remains approximately 
constant, but the potential energy of the gas decreases with 
p* , and so the escape of the gas is favored. 

The effect of rotation is shown in Fig. 2, again for 
the isotropic rotator and for the maximum rotation solution 
(with 7 = 1), for various flattening of the luminous matter. 
The variation in % is slightly larger than in the case of no 
dark matter, especially when the dark halo is spherical. In 
any case, this variation for xn — 1 is always very small, and 
so the result of the case of no dark matter is confirmed: the 
effect of rotation is negligible in the energy balance of hot 
gas flows. 

6 DISCUSSION AND CONCLUSIONS 

In this paper we investigate which flow phase of the hot 
gas is most likely in stellar systems characterized by various 
degrees of flattening and rotation, for the WOI scenario. 
Energetic arguments have been used that involve the ratio 
X between the energy required to extract the hot gas from 
the galaxy potential well and the energy supplied for this 
process. For a general stellar system, we proved analytically 
that the critical variations in % (i-e., from > 1 to < 1 or 
viceversa) can be produced only by a change in the galaxy 
structure; rotation - if not fully thermalized - acts in the 
sense of increasing a % > 1 and decreasing a % < 1 . 

Realistic galaxy models for SOs - with and without dark 
matter - have been built by using the Miyamoto-Nagai den- 
sity distribution, and by comparing their central stellar ve- 
locity dispersion (jp with that of a spherical galaxy of the 
same ig . The analytical formulae for the virial properties 
of the two-component galaxy models have been derived, for 
the Satoh's ^-decomposition and for the maximum rotation 
solution, together with the projected edge-on ordered and 



random stellar velocity field, in the case of no dark matter. 

It is shown that reasonable flattenings of the galaxy 
structure at fixed ig can change a % > 1 of a round system 
into a X < 1, i.e., can make the gas less bound (for a (jp of 
the flat system comparable to that of the round one). The 
role of rotation for the adopted galaxy models has proved 
to be negligible, because it changes % by only a few percent 
even in the maximum rotation case. So, the problem of the 
unknown amount of thermalization of the ordered motions 
does not affect our conclusions. 

The efficiency of the fiattening and the inefficiency of 
rotation in changing the nature of the gas fiow phase is con- 
sistent whith the observational result that - at fixed ig - 
ix/^B is on average lower not only in SO galaxies with re- 
spect to Es, but also in flat ellipticals with respect to round 
ones; evidence of this trend has in fact been found by Es- 
kridge et al. (1995b). 

Other possibilities could be imagined to explain the X- 
ray underluminosity of SOs with respect to Es. One could be 
a higher SNIa rate, due to the possible presence of a younger 
stellar population in the disk: this could favor the escape of 
the gas. This solution does not seem likely, in view of the 
fact that also fiat Es show the same trend, and it cannot 
be sustained that the stellar population of fiatter Es has 
different properties from that of rounder Es. Similarly, the 
possible presence of an X-ray absorbing cold material in the 
disk of SOs does not seem to be a general explanation of the 
effect . 

We conclude pointing out a problem with the purely 
energetic approach of the present investigation: x is a glob- 
al parameter, and so it does not give much indication when 
there is a sort of "decoupling" in the fiow between different 
regions of the galaxy. This is the case, for example, of the 
partial wind phase (see, e.g., Pellegrini &; Fabbiano 1994). 
In connection with this point we also remark that the re- 
liability of X has been proved by CDPR just in the case 
of spherical models, and the above mentioned "decoupling" 
could be more likely for fiattened systems. Hydrodynam- 
ical simulations are required for a definite answer, and an 
extension of the CDPR study to axisymmetric models is in 
progress. 
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APPENDIX A: PROJECTED PROPERTIES 
OF MN MODELS 

In this Appendix the analytical expressions for the projected 
velocity dispersion associated with a stellar model described 
by the density-potential pair (9)-(10) are given. The pro- 
jection of the velocity dispersion is performed both for the 
Satoh and the maximum rotation decomposition. 

We start evaluating the projected surface density. 
Defining 



I arcos (s)/Vl - s2 ifO<s<l; 
F{s) = \l i{s = l; 



(Al) 



L arcch(s)/-\/s2 — 1 if s > 1 

the face-on projection in the center is given by Sh.(0) 
{M*/2Tvb'^)E*{s), where 



~ 2 + s^-3sF{s) 

^H")- (1_,2)2 ' 



(^2) 



with Sh.(1) = 2/5. EH.(ii) could not be determined analyti- 
cally outside the center. 

The edge-on projection of the density at the point (ii, z) 
in the projection plane is given by 

Mn-fe^ aR^ + {a + 2C){a + Cf 



2-K C3[ii2 + (a + C)2 



(^3) 



where C = Vz^T^ ■ This expression was already obtained 
by Satoh (1980), and is given here for completeness. 

Moving to the projection of the squared velocity, this is 
written as 



-C 



p^< v,n >2 dl, 



where <,> is the inner product, n is the l.o.s. direction, I 
is a parameterization of the l.o.s. integration path, and the 



overline represents the mean over the phase space (see, e.g., 
BT87, p. 208). Expanding (A4) one obtains = Vp + a^, 
where Vp is the projection of the square of the component 
along n of the ordered velocity, and (jp is the projected ve- 
locity dispersion. Note that (jp is not the line-of-sight (i.e., 
the observed) velocity dispersion if a net rotation is present; 
this is instead given by: 

/oo 
p*{< v,n > -■wios)^ dl = Eh.(up - uios), (^5) 
- oo 



where 



f 

J — C 



/?=(=< v,n > dl. 



{A6) 



The analytic solution of the integration in (A6) could not 
be found. Note however that for the central face-on and 
edge-on projections it results a^^^ = (jp. 

We now evaluate the two components of Vp = Vp + 
(Jp for an edge-on projection. The integration in (A4) is 
made by changing the angular coordinates of n to the radial 
coordinate; in the case of the Satoh decomposition, using 
(13) one obtains 

V^{R,z) = k''vUR,^), {A7) 



and 

al{R,z) = S?{R,z) + {l-k^)vUR,^), 



{A8) 



where the meaning of the two new functions 5jg and Vjg is 
evident. After integration 



GM'iah'^iR^ 



32C3[ii2 + (a + C)2]5/2' 



and 



2 _ GM'ih'^i{a + Cf 

64C2[ii2 + (a + C)2]5/2' 



(^9) 



Uio) 



In the case of the maximum rotation decomposition, apply- 
ing (15) one obtains: 



V^{R,z) = yUR,z) + Sl{R,z), 
and 

a-l{R,z) = Ss{R,z), 

with 5c + 5s = 5;g. After integration. 



(An) 

{A12) 



S*5c 



r 8R' + 20{a + C)'R^ + 15{a + C)''R \ 
^ |« [ii2+(, + ^)2]5/2 I ' (^13) 



and 



2 GMSb^ 



64C2(a + C)4 



j 8R'' + 12{a + C)'R'+3{a + C)'' \ 
1 [il2+(a + C)2]3/2 

Note that, because of the intrinsic underdeterminacy of the 
Jeans equations, Vp has to be independent of the particular 
adopted velocity decomposition, and a check shows that in 
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both of our cases Vp = V[g + 5;^. 

Finally we compute the face-on projection of the veloc- 
ity dispersion: after integration along the z axis one obtains 
EH.(0)cr2(0) = {GMS /2Tvb^)E*{s)al{s), where 

Eh.(s)?p(s) = ^ + 



3(8^'' - Ss'^ + 7s'' - 2)F{s) - s{26s'^ - 17s'' + 6) 



, (^15) 



12s4(l -s2)3 

with Eh.(0)?2(o) = 37r/32, and Eh.(1)?p(1) = 7r/4 - 16/21. 

APPENDIX B: VIRIAL THEOREM FOR 
MN+MN MODELS 

We compute here the analytical expressions for the terms en- 
tering the virial theorem of a MN density distribution that 
is interacting with another MN distribution. In order to 
simplify the notation and the resulting formulae, all the en- 
ergies appearing in the virial theorem are normalized to the 
reference energy £ = GMS/8b. From the linearity of the 
Jeans equations with respect to the potential, each of the 
contributions to the virial theorem is split in a part involv- 
ing the interaction of one MN distribution with itself, and 
another involving the interaction with the other MN compo- 
nent; so \W\ = |VFH.| + |VFh|,T = TH.+Th andn = nH.+nh. 
We start studying \W\, that is obviously independent of the 
assumed velocity decomposition: 



i/ 



p*^*dV = £W*{s), 



and after integration over the whole space, one gets 



W^s) 



1 - 2s'' 



s(l-s2) 2s2 s2(i_52) 



(Bl) 



(52) 



where J'(s) is given in (Al), and so VFh.(0) = 37r/4, and 
VFh.(1) = 8/3 — 7r/2. The interaction energy can be written 
as 



), (S3) 



with 



Wy,{c,s) 



\ (1 - c^Y (2 5\ ^, , 



; (54) 



so W^h(0,0)/7J = 37r/4, and W^h(l,0 < s < 2)/7J = 7r/2 - 
16/15 + s(56/15 - tt). 

The gravitational energy of the stellar distribution in- 
side the dark matter halo potential is given by: 



i[/hi=- y p*$h<i^"=ft/h(c), 

which after integration becomes 
[7h(c) = 27JW^H.(c) 



(55) 
(56) 



where the function on the r.h.s. is the same given in (B2). 
Moving now to the kinetic virial terms, and defining 



T = 
and 



1 I 



dV, 



His = J p*cr^dV, 



(57) 



(58) 



after a simple algebraic manipulation of (13) and (14), we 
obtain for the Satoh's decomposition: 



T = £k''T, 
and 

n = f[3n;, + 2(i-fc2)f]. 

For the maximum rotation, 

T = f(f + 0.55;,), 

and 

n = £2liu. 



(59) 
(510) 
(511) 
(512) 



So the analytic expressions for the two functions T and Xlig 
are needed. Integrating by parts on z, we have 
/»oo /»oo 

Hi, = 47r y zdz J Rp* -^'^^ = His,* + n;,,h, {BU) 

and for the ordered kinetic energy, integrating by parts on 
R and then on z: 



T = 2 



/»00 PC 

■K I dz I 

Jo Jo 



R^P* ^^dR - Hi, = To - n;,. (514) 



After tedious calculations, one obtains: 



nis,H.(s) 



TT 1 (2a^ - l)F{s) 

2I2 " s(l - s2) + s2(l-s2) 



(515) 



with nis^H.(0) = 7r/4, and nis^H.(l) = i"/2 - 4/3. A formally 
identical integration gives: 



nis,h(c) = 7Jnis^H.(c). 

Finally, for the ordered component Tq, we have 

~ 1 TT F{S) 

with fo,H.(0) = 7r/4, fo,H.(l) = 2 - 7r/2, and 
(i^fo,H(c,.) = -f +(cVi)i^(c)+ 



(516) 
(517) 



2\2 



7r(l -c^) 



,(518) 



with To,h(0,0)/7J = 7r/4, and To,h(l,0 < s < 2)/7J 
2/15 + s(28/15 - 7r/2). 



